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Several model fluids in contact with planar, spherical, and cylindrical walls are investigated for 
small number densities within density functional theory. The dependence of the solid-fluid interfa¬ 
cial tension on the curvature of spherical and cylindrical walls is examined and compared with the 
corresponding expression derived withiir the framework of morphometric thermodynamics. Partic¬ 
ular attention is paid to the implicatioirs of the choice of the interface location, which underlies the 
definition of the interfacial tension. We And that morphometric thermodynamics is never exact for 
the considered systems and that its quality as an approximation depends sensitively on the choice 
of the interface location. 
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I. INTRODUCTION 

In recent years morphometric thermodynamics has 
been used to understand the influence of geometric con¬ 
straints on the thermodynamic properties of fluids [1-5]. 
This approach has been motivated by a theorem of inte¬ 
gral geometry which is often referred to as “Hadwiger’s 
theorem” [6, 7]. For example, in Ref. [1] it has been ap¬ 
plied to a fluid of hard spheres bounded by a hard wall. 
Based on the assumption of being an additive, motion- 
invariant, and continuous functional of the shape of the 
walls, the grand canonical potential Deq of a confined 
fluid is given, in accordance with Hadwiger’s theorem, by 
a linear combination of only four geometrical measures 
which characterize the shape S of the confining walls: 
volume V, surface area A, integrated mean curvature C, 
and Euler characteristic Y: 

Deq[^] = -pV[S] + 7oM[5] + kC[S] + RY[S]. (1) 

According to Eq. (1) the pressure p, the interfacial ten¬ 
sion for a planar wall 70, and the bending rigidities k 
and K, unlike the geometrical measures, do not depend 
on the shape of the bounding container. This structure 
turns morphometric thermodynamics into a very attrac¬ 
tive tool, because once the thermodynamic coefficients 
p, 7 o,K, and R have been determined — preferably by 
considering simple geometries S — the grand canoni¬ 
cal potential can be readily calculated even for systems 
bounded by complex geometries the shape S of which en¬ 
ters Eq. (1) only via the measures V,A,C, and Y. This 
strategy has been used, e.g., in Refs. [2, 3, 5] in order to 
calculate solvation free energies of proteins with complex 
shapes. 

As a thermodynamic quantity, which follows from the 
grand canonical potential, the interfacial tension 7 ac¬ 
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quires the simple morphometric form [1] 

7 = 7q -I- kH -I- RK, (2) 

where H = C/A and K = Y/A denote the averaged 
mean and Gaussian curvatures, respectively, of the con¬ 
fining wall [1]. In the case of geometries with constant 
curvatures, i.e., spherical and cylindrical walls with radii 
of curvature R, Eq. (2) can be written as 

70 + ^ + ^, spherical wall 
7o-l-^, cylindrical wall, 

where 71,2 denote coefficients independent of the curva¬ 
ture 1/R and where 72 = 0 in the case of cylindrical 
walls. 

On one hand morphometric thermodynamics has al¬ 
ready been applied to several physical systems confined 
to geometries with complex shapes [2, 3, 5] and it has 
been found to be in agreement with certain data of hard 
spheres [1, 4, 8] and of hard rods [9] in contact with hard 
walls. On the other hand, recently evidences have been 
provided that a description as in Eq. (2), where the cur¬ 
vature expansion of the interfacial tension for a spherical 
wall terminates after the quadratic and that of a cylindri¬ 
cal wall after the linear order in 1 /R, could be incomplete 
[ 10 - 12 ]. 

Here we propose an explanation for the observation 
that for certain studies morphometric thermodynamics 
appears to be applicable, whereas it is not for others. To 
that end, we analyze several model fluids with small num¬ 
ber densities in contact with curved walls using density 
functional theory (DET) within the second virial approx¬ 
imation (see Sec. H). This technically simple approach 
allows us to study various kinds of interactions among 
the fluid particles (Eq. (4)-(6)) with high numerical pre¬ 
cision for low densities (see, e.g., Eig. 3). For this reason 
we abstain from considering fluids with high densities as 
they occur, e.g., at two-phase coexistence. It turns out 
that those basic features of morphometric thermodynam¬ 
ics we are aiming at reveal themselves already at low den¬ 
sities for which one can control the numerical accuracy 
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sufficiently in order to address reliably the correspond¬ 
ing issues. Moreover we present exact results for an ideal 
gas fluid confined by non-hard spherical and cylindrical 
walls. We focus on the interfacial tension 7 and compare 
its curvature dependence with the one predicted from 
morphometric thermodynamics (Eq. (3), see Sec. III). It 
turns out that the morphometric form of the interfacial 
tension is indeed not valid exactly and that its quality as 
an approximation depends on the choice of the location 
of the interface underlying the definition of the interfacial 
tension (see Sec. IV). 


II. MODEL 

We consider a simple fluid composed of particles which 
interact via an isotropic pair potential U (r) which is char¬ 
acterized by an energy scale Uq, a length scale L, and, 
for technical convenience, a cut-off length Lc such that 
U{r > Lc) = 0. We focus on three distinct types of pair 
potentials: 


The geometrical properties of the walls enter the de¬ 
scription only via the external potential V®’'*. In order to 
determine the thermodynamic coefficients in Eq. (I) we 
consider the bulk fluid and the fluid in contact with pla¬ 
nar, spherical, and cylindrical walls exhibiting constant 
curvatures. For these choices the densities in Eq. (7) de¬ 
pend on a single spatial variable only. 

• The uniform bulk corresponds to a spatially con¬ 
stant external potential, which can be set to zero 
without loss of generality: 

l3V^^\r)=0. (9) 

As a consequence the equilibrium density 0 eq(r) = 
is independent of the position r. 

• A planar wall leads to an external potential 


j3V‘^^\v) 


00 , z < 0 
/3V(z), z>0. 


( 10 ) 


• the square-well (C/q < 0 ) or square-shoulder (I/q > 
0 ) potential 

pU{r < Lc) = m (4) 

with L = Lc, 

• the Yukawa potential {Uq > 0) 

f3Lf{r < Lc) = l3Lfo^exp(^-^'j , (5) 

• the Lennard-Jones potential {Uq > 0) 

mr<Lc)=PUo(^(^^y (6) 

In the following non-uniform number density profiles 
p(r) of the fluid in contact with walls are studied. To 
this end density functional theory [13] offers a particu¬ 
larly useful approach. Since the present investigation is 
focused on low densities, we use the density functional 
n[(?] within the second-virial approximation [14, 15]: 

j3Vl[Q] = Jd^r g{r) (in (p(r)T^) - 1 - /3/i -|- /3V®’'‘(r)) 

( 7 ) 

+ ^ dV' ei(r) 0 (r') (1 - exp(-/?C/(|r - r'|))), 

where /3 = 1 /(A:bT) is the inverse temperature, V®’'*(r) is 
the external potential due to walls, and p, is the chemical 
potential p. shifted by a constant (with respect to r and 
g) given by the thermal de Broglie wavelength A and 
the length scale L of the pair potential among the fluid 
particles: 



This potential implies that the equilibrium density 
geq{z) is identically zero for z < 0. Therefore it is 
not possible for the centers of the fluid particles to 
get closer to the wall than z = 0. In the following 
we call the set of accessible points of the centers 
of the fluid particles closest to the geometrical wall 
surface the reference surface. In Eq. (10) V{z) rep¬ 
resents the excess part of the wall potential (i.e., in 
excess of V®’'‘(z < 0 ) = 00 ) as a function of the dis¬ 
tance z from the reference surface. Depending on 
the wall-fluid interaction potential, the geometrical 
wall surface at position Xg in Fig. I and the refer¬ 
ence surface at position X in Fig. 1 can be distinct. 

• A spherical wall with radius R of the reference sur¬ 
face is characterized by an external potential 


/3 V®’'* (r) 


00 , r < R 
PV{r), r > R, 


( 11 ) 


where V(r) represents the excess part of the exter¬ 
nal potential as a function of the distance r to the 
center of the sphere. 

• A cylindrical wall with radius R of the reference 
surface is characterized by an external potential of 
the same form as the one in Eq. (II); however, 
in this case, r and R measure distances from the 
symmetry axis of the cylinder. 

In accordance with the variational principle underly¬ 
ing density functional theory [13] the equilibrium density 
Peq(r) minimizes the functional in Eq. (7). The corre¬ 
sponding Euler Lagrange equation 

5 pm 


Pg := Pp — 3 In 


( 12 ) 
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FIG. 1: The reference surface (solid line) at position X is 
determined by the external potential. It consists of the set 
of all reachable locations of the centers of the fluid particles 
which are closest to the geometrical wall surface (dash-dotted 
line) at position Xg. Depending on the wall-fluid interaction 
potential, the geometrical wall surface and the reference sur¬ 
face can be distinct. The interface of area A (dashed line) — 
with respect to which the interfacial tension 7 in Eq. (16) is 
defined — may differ from the reference surface. The param¬ 
eter 5 measures the offset of the interface with respect to the 
reference surface. This sketch refers to the cases of spherical 
or cylindrical walls for which X = R. However, the concept 
involving the parameter 5 is valid also for other geometries, 
in particular for a planar wall with X = z = 0. If the dashed 
line runs within the interior of the reference surface, 5 is taken 
to be negative. The fluid volume V, which refers to the set 
of all points being not closer to the geometrical wall surface 
than the interface, is shaded in grey. 


is solved numerically. From the equilibrium density pro¬ 
file Peq(r), the grand canonical potential fleq follows from 
[13] (Eq. (7)) 


fleq = fulfil = Peq]- (13) 

We have verified that the hard wall sum rule (see, e.g., 
Ref. [ 8 ]) is fulfilled by the functional in Eq. (7). 


III. DISCUSSION 


The interfacial tension 7 is defined by 


/lOeq - 

AL-^ 


/3fteq + PpV 
AL-^ 


(14) 


as the work Oeq — fleq"' interfacial area A required 
to create the interface [16]. On the right hand side of 
Eq. (14) the relation 

= _Pp V (15) 


between the grand canonical potential Ogq"', the fluid 
volume V, and the bulk pressure p has been used [16]. 

The wall-fluid interfacial tension 7 is not an observable 
because according to Eq. (14) its value depends — via V 
and A — on the arbitrary choice of an interface position. 
In order to characterize various interface conventions we 
introduce a parameter S which measures the offset of a 
chosen interface position with respect to the reference 
surface (see Fig. 1). In addition to the interfacial area A 
the choice of an interface position also determines what 
is called the fluid volume V, which refers to the set of 
all points being not closer to the geometrical wall sur¬ 
face than the interface. As a consequence A and V are 
functions of X -|- 5 where X characterizes the reference 
surface (see Fig. 1). On the other hand Oeq depends on 
X only, because due to Eq. (13) only the parameters of 
the substrate potential (i.e., X) enter Oeq. Accordingly, 
one has 




/30eq(A) + Pp V{X + 6) 

A(X -b S)L-^ 


(16) 


Our data are calculated within the convention 6/L = 0 
as this choice is convenient for various types of interac¬ 
tions U. Moreover, in certain studies (see, e.g., Ref. [17]) 
this choice has been argued to be “the natural one from 
the point of view of statistical mechanics”. Eigure 2 
shows the interfacial tension 7 for hard spheres in con¬ 
tact with spherical (Fig. 2(a)) and cylindrical (Fig. 2(b)) 
walls with various radii R of the reference surface. The 
packing fraction 77 = Trg^^^^L^/6 is chosen sufficiently 
small such that the second virial approximation is valid. 
The plot for the cylindrical walls (Fig. 2(b)) reveals a 
non-linear increase, similar to the case of spherical walls 
(Fig. 2(a)). This means that the curvature expansion of 
7 in terms of powers of 1 /i? does not terminate after the 
hrst-order term, in contradiction to the prediction from 
morphometric thermodynamics for cylindrical walls. 

Evaluating Eq. (16) for (5 = 0 and arbitrary (5 yf 0 and 
exploiting that fleq does not depend on 6 leads to 


l2/57(A,5) S = 0) 


(17) 


According to Eq. (17) the interfacial tension 7 (A, <5 = 0) 
calculated for the convention d/L = 0 can be translated 
to that for any other choice of the convention. For ex¬ 
ample the convention 6/L = —0.5 is often used when 
discussing hard spheres confined by hard walls because 
this choice renders the interface to coincide with the ge¬ 
ometrical wall surface at position Xg in Fig. 1, which, 
in this case, is separated from the reference surface at 
position X in Fig. 1 by a distance given by the particle 
radius L/ 2 . 

In Fig. 3 simulation results of Laird et al. [4] for 
hard spheres with packing fraction 77 = 0.02656, obtained 
for 5/L = —0.5, are plotted together with the data of 
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FIG. 2: The interfacial tension 7 of hard spheres with di¬ 
ameter L for the convention (5/L = 0 as a function of the 
radius R of the reference surface forming the boundary of a 
spherical (a) and a cylindrical (b) hard body. Polynomials 
(lines) as predicted by morphometric thermodynamics (up to 
order 1/{R + SY for spherical walls and up to order l/{R + 5) 
for cylindrical walls) have been fitted to the numerical data 
(dots). In panel (a) the fit shows good agreement with the 
data, whereas panel (b) clearly shows deviations from the 
predicted behavior. This means that for the cylindrical wall 
the curvature expansion of 7 Eq. (3) does not terminate after 
the first-order term, in contradiction to the prediction from 
morphometric thermodynamics. The hard spheres interact 
among each other via the square-shoulder potential U{r) in 
Eq. (4) with pUo = 00 . The chemical potential is chosen as 
Pfi. = —2.768839 which corresponds to the low packing frac¬ 
tion rj = f « 0.02656. 


Fig. 2 which have been translated into the convention 
6/L = —0.5 according to Eq. (17). The agreement with 
the simulation data is very good. In particular, within 
this convention for 6, in the case of a cylindrical wall 
(Fig. 3(b)) the data points almost coincide with a straight 
line, in accordance with the prediction of morphometric 
thermodynamics. This finding has also been confirmed 
for the packing fractions ry Ri 0.053 and O.IOI, for which 
the respective plots are qualitatively similar to those in 


FIG. 3: Comparison of the simulation results of Laird et al. [4] 
(squares □) with the corresponding data of the present DFT 
study (crosses x) for the packing fraction 7 ss 0.02656. The 
latter have been obtained by translating the data displayed 
in Fig. 2 (for 5 = 0) into the convention 5/L = —0.5 by using 
Eq. (17). In contrast to Fig. 2(b), the plot of the interfacial 
tension in the case of a cylindrical wall (b) almost coincides 
with a straight line. 


Fig. 3 except that, as expected, the results of the present 
virial expansion deviate more and more from the simula¬ 
tion data upon increasing the density. 

Figures 2 and 3, which are based on the same micro¬ 
scopic system, show that the interfacial tension depends 
strongly on the choice of the convention for <5. Upon 
varying 5 not only the sign of 7 may change, as already 
noted in Ref. [ 8 ], but also the magnitude and even the 
qualitative functional form, which is revealed clearly in 
the case of cylindrical walls. 

Figure 4 shows the behavior of the same microscopic 
system as above (hard spheres exposed to a hard cylin¬ 
drical wall) described in terms of two conventions for 
T] « O.OI. The data are presented in log-log plots which 
facilitate the identihcation of power laws in I/(i? -|- 5). 
In this presentation the contribution 70 (5) of the planar 
wall is subtracted so that the plotted quantity vanishes 
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FIG. 4: Dependence of the interfacial tension 7 for hard 
spheres at a cylindrical hard wall on the radius R of the 
reference surface and on the shift parameter S. The con¬ 
tribution of the planar wall 70 (5), which corresponds to the 
limit i? —>■ 00 , has been subtracted. The numerical data (red 
dots) are obtained within the convention 5/L = 0 (a) and 
translated to the convention S/L = —0.5 (b) via Eq. (17). 
The three straight lines represent the terms n = 1,2,3 of 
the curvature expansion in Eq. (18). In the case 5/L = Q 
the coefficient 72 ( 5 ) is the largest one, i.e., for small R the 
contribution j3'y2{5)L*/[R + 5)^ is the dominant one so that 
a deviation from the leading behavior ~ l/(i? + 5) becomes 
obvious for R/L < 3. Within the convention 5/L = —0.5 
the leading coefficient 7i(5) is the largest one, i.e., the lead¬ 
ing term /3y-i{5)L^/{R + 5) is the dominant one. Here U(r) 
is given by Eq. (4) with [/q —>■ 00 and = —3.88 so that 
1 ] « 0.01. Note that at {R-\- 5)/L = 1 the plots render the 
values of the dimensionless coefficients /3'yn{5)L‘^. 


for R —> 00 . Within the convention S/L = 0, at R ^ 3L 
there is a crossover between two power laws. Thus the 
dependence of the interfacial tension on 1 /i? consists of 


more than the leading term ^ 1 /i? which, according to 
morphometric thermodynamics, would be the only one 
allowed for the cylindrical configuration. The behavior 
is different within the convention S/L = —0.5. There, 
within this presentation, the interfacial tension is repre¬ 
sented by an almost straight line throughout the whole 
range of {R + S)/L shown. In order to analyze the curva¬ 
ture dependence of the interfacial tension more quantita¬ 
tively, we assume that 7 ( 1 ?, (5) can be expanded in terms 
of a power series in l/(i? -I- <5): 


L^j3-i{R,S) = L^ 




n—Q 


liln{S) 

{R + sy 


(18) 


For various radii R of the reference surface of the curved 
wall, the interfacial tension 7 is calculated numerically 
(red dots in Fig. 4) and fitted to the curvature expansion 
Eq. (18) with n < 10. This way the coefficients 'Yn(S) 
have been determined. Only the coefficient 70 (i5), which 
is the interfacial tension at a planar wall, can be ob¬ 
tained independently without fitting. In Fig. 4 the terms 
L^~^^yn(S)/(R + S)^ corresponding to n € {1,2,3} in 
the curvature expansion of Eq. (18) are plotted as lines. 
Within the convention S/L = 0 the quadratic coeffi¬ 
cient 72 , which vanishes within morphometric thermo¬ 
dynamics, is even larger than the linear coefficient 71 ; 
this explains the crossover between two power laws de¬ 
scribing the numerical data (red dots). However, within 
the convention S/L = —0.5 the first-order coefficient 71 is 
much larger than the higher order coefficients 7 „,n > 2 , 
so that in this case morphometric thermodynamics is 
a very good approximation of the exact curvature de¬ 
pendence of the interfacial tension. Within the class 
of systems with square-well-like or square-shoulder-like 
particle-particle interactions, we have studied a large va¬ 
riety of configurations within the convention S/L = 0 in 
the same way as shown in Fig. 4. Fluids with packing 
fractions ry S {0.01, 0.02, 0.05, 0.10} and with interac¬ 
tion strengths pUo € {— 0 . 1 , 0 . 1 , 1 , 00 } have been exam¬ 
ined near spherical and cylindrical hard walls. In addi¬ 
tion to the interfacial tension 7 , the dimensionless excess 
adsorption F [16] 


L^iX.S) 


N{X)-0y^V(X + S) 

A{X + S)L-^ 


(19) 


has been calculated for Uq > 0 where N{X) = 
fv(x) denotes the number of fluid particles. Our 

corresponding observations can be summarized as fol¬ 
lows: 


• Apart from opposite signs, both the excess adsorp¬ 
tion F and the interfacial tension 7 exhibit a similar 
dependence on the radius of curvature R. 

• The functional form of (/ 37 (i?, 0) —/ 37 o( 0 )) is 
similar when comparing fluids with the same bulk 
state near a spherical and a cylindrical hard wall. 
Because in all considered cases the third order coef¬ 
ficient 73 is smaller than 71 and 72 , morphometric 
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thermodynamics turns out to be a better approxi¬ 
mation for spherical than for cylindrical walls. 

• Decreasing packing fractions r] or interaction 

strengths C/q > 0 result in a shift of the crossover 
between the power laws describing — 

7o(0)) towards larger radii R, i.e., the second coef¬ 
ficient 72 becomes larger in comparison to the first 
coefficient 71. Apart from that, the behavior of 

{(3^{R,0) — /?7o(0)) is similar to the one shown 
in Fig. 4. For Uq < 0, L'^ {(3^{R,0) — /?7o(0)) ex¬ 
hibits a zero because the coefficients 71 and 72 have 
opposite signs. 

• For some of the above systems the data have been 
translated to the convention S/L = —0.5 by us¬ 
ing Eq. (17). In these cases, for the spherical 
wall configurations the third order coefficient 73 is 
much smaller than 71 and 72 , and for the cylin¬ 
drical wall configurations the leading coefficient 71 
is much larger than the subleading ones. There¬ 
fore, within this convention morphometric thermo¬ 
dynamics turns out to be a very good approxima¬ 
tion. 

The comparison of the plots in Figs. 4(a) and (b) shows 
that the coefficients 7 n(iJ) in the curvature expansion 
(Eq. (18)) indeed depend on the chosen convention 6. 
In order to examine the implications of shifting the po¬ 
sition of the interface (see Fig. I) this dependence is in¬ 
vestigated more closely. The derivative of Eq. (16) with 
respect to 6 for fixed R leads to 

L^L2^7(i?, 5) = -L^Pp - dj^L^PjiR, S) 

^ (2, spherical wall 

I 1 , cylindrical wall. 


where V'{R + 6) = —A{R + 6) and A'{R + 6)/A{R + S) = 
d/{R -I- 6) have been used. The derivative of Eq. (18) 
with respect to S gives 



L^Pl{R, S) = 


00 


E 

n—O 


f.n 

1 {R + 5Y 


- L^n 


PlnjS) \ 

(i? + 5)"+i/■ 


( 21 ) 


Equating Eqs. (20) and (21) and using Eq. (18) leads to 


-l^Pp = 

X 


l^Pi'o{S) + Y1 

n—1 

{L^P^'^{6) + {d 


Jr+W 

-n + 1 )L^/ 37 „_i(( 5 )} 


( 22 ) 


for all R. Comparison order by order in {R + S) ^ in 
Eq. (22) renders 

0((i? + <5)°): -L^Pp = L^pj'oiS) (23) 


and 

0{{R + 6)--),n>l: 

L^PlUS)+{d-n + l)L'^P^n-i{S) = 0. 

Integration of Eqs. (23) and (24) with respect to <5 leads 
to the following iterative algorithm for determining the 
dependence of the coefficients 7 „( 5 ),n > 0 , on (5: 


n = 0 : 

= L^Ploi^) - 

n > 1 : 


s _ 

= T^/ 37 „( 0 ) + {n-d-l) J ^T^/ 37 „_i(( 5 ) 

0 



spherical wall 
cylindrical wall. 


(25) 


The dependence of the coefficients 7 n(^) on S is fully 
determined once in Eq. (25) the values of 7 n'( 0 ) for all 
n' <n are known. Here the values 7n(0) are obtained by 
fitting Eq. (18) for n < 10 to the numerical data within 
the convention 6/L = 0. Considering terms in Eq. (18) 
to such high orders was necessary in order to achieve a 
sufficiently high precision for the actually interesting co¬ 
efficients 7 i(( 5 ), ..., 73 (^) (see Eigs. 4-9); taking the addi¬ 
tional terms of order n > 4 into account guarantees that 
these coefficients 71 (<5),..., 73 ( 5 ) are not affected by the 
fast-decaying contributions of the full curvature expan¬ 
sion. Thereby we have found that the ratio of the leading 
coefficients for the spherical wall, 7 f ( 0 ), and for the cylin¬ 
drical wall, 7 i( 0 ), takes the value 7 i( 0 )/ 7 f( 0 ) = 2 for all 
systems considered here, with a relative deviation of 10 “^ 
or less. Comparison of that numerical result with the ex¬ 
act relation 7i(^)/7i(<5) = 2 , which follows from the fact 
that the total curvatures J = 1/Ri + l/i ?2 of a sphere, 
J®, and of a cylinder, J'^, are related by J® = 2J‘^ (see 
Refs. [I, 10]), validates our numerical approach. Mo¬ 
tivated by this relationship between the leading coeffi¬ 
cients we have considered also the ratio of the next-to- 
leading coefficients. Eor small packing fractions we have 
obtained 72(<5)/72('5) ~ 8/3 for <5 = 0, independently of 
the particle-particle interaction potential U. The value 
of this ratio does depend on the convention for <5 because 
the spherical coefficient 7 |(( 5 ) varies with 5, whereas the 
cylindrical coefficient yf is constant (see discussion be¬ 
low). Actually, the value 8/3 can be obtained from the 
exact analytical expressions for the surface tension in the 
low-density limit for a fluid of hard spheres [11]. More¬ 
over, the exact expression in Eq. (A7) (see Appendix A) 
describes the deviation of the ratio 72(0)/72(0) from 8/3 
for an ideal gas of particles as function of the strength 
PV^"^ of a short-ranged excess external potential in ad¬ 
dition to the hard wall potential. 

It is interesting to pay special attention to the coeffi¬ 
cient ^n=d+i in Eq. (25) which is the coefficient of the 
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FIG. 5: Reduced coefficients 'yn{5),n £ {0,1,2, 3}, charac¬ 
terizing the curvature expansion in Eq. (18) as function of 
the shift parameter 5 and as obtained from Eq. (25). The 
data correspond to a hard sphere fluid at /3/r = —3.88, so 
that T] « 0.01, exposed to hard spherical (a) or cylindrical (b) 
walls. 


lowest order being not in accordance with morphomet¬ 
ric thermodynamics. For this order n = d 1 one has 
n — d — 1 = 0 and therefore 7 d+i is constant in 5 (see 
Eq. (25) for n > 1). This checks with Fig. 4 (correspond¬ 
ing to d = 1 ), where the values of the coefficients 72 can 
be read from the lines 1/(7? + d)^ at {R + 5)/L = 1. 
The value of 72 is not vanishing and it is the same in 
both conventions for S. This implies that within mor¬ 
phometric thermodynamics the curvature expansion is 
not exact for any convention for S. On the other hand, 
if, as a consequence of Eq. (25), the i?-dependence of 
7 ( 7 ?, S) within morphometric thermodynamics would be 
exact for any single convention for S, it would be exact 
for all conventions for S. However, this statement is of no 
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FIG. 6: Reduced coefficients 'yn{S),n G {0,1, 2, 3}, charac¬ 
terizing the curvature expansion in Eq. (18) as function of 
the shift parameter S and as obtained from Eq. (25). The 
data correspond to a fluid with a square-shoulder pair poten¬ 
tial U{r) (Eq. (4)) with /3Uo = 0.1 and /3/i = —3.95 so that 
rj ~ 0.01. The fluid is exposed to hard spheres (a) or hard 
cylinders (b). 


practical use, because, on the basis of numerical data, it 
is virtually impossible to prove that there is a convention 
for d within which the morphometric form of the interfa¬ 
cial tension is exact. In contrast. Fig. 4 shows that even 
if the non-morphometric coefficients jn(d) (n > 2 for a 
cylindrical wall) are numerically small for one convention 
for S (see Fig. 4(b)) they may be large for another (see 
Fig. 4(a)). The reason for this is that the operation of 
approximating the curvature-dependence of the interfa¬ 
cial tension by the form predicted within morphometric 
thermodynamics does not commute with the operation 
of shifting the interface. 

Figure 5 shows the dependence of the coefficients 'yn(d), 
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FIG. 7: Same as Fig. 6 for a fluid with a Yukawa pair potential 
(Eq. (5)) with pUo = 0.1, /3/x = —3.94, and Lc/L = 5 so that 

ri « 0 . 01 . 



S/L 



SjL 

FIG. 8: Same as Fig. 7 for a fluid with a Lennard-Jones pair 
potential (Eq. (6)) with fiUo = 0.1, /3/i = —3.92, and Lc/L — 
5 so that rj « 0.01. 


0 < n < 3, of the curvature expansion Eq. (18) on the 
shift parameter <5. Except for the first non-morphometric 
coefficient (i.e., 73 for a sphere and 72 for a cylinder), 
which is constant in 6, the coefficients 7 ™ (<5) vary over 
several orders of magnitude upon changing S. This is 
particularly pronounced near S/L = 0, where the mor¬ 
phometrically allowed coefficients are small; in the case 
of a cylindrical wall 71 is even smaller than the leading 
non-morphometric coefficient 72 . However, apart from 
this region around S/L = 0, e.g., at S/L = —0.5, the 
morphometrically allowed coefficients exceed the leading 
non-morphometric coefficient by several orders of mag¬ 
nitude. These observations are in agreement with the 
findings of Fig. 4 which is based on the same system and 
where for each convention for S the coefficients, the values 
of which are rendered a.t {R +S)/L = 1, have been fitted 


independently. Figure 5 demonstrates that the interfacial 
tension cannot be represented exactly by the form ob¬ 
tained within morphometric thermodynamics, and that 
the quality of the approximation of the interfacial ten¬ 
sion by the morphometric form depends on the position 
of the interface parameterized by the shift S. 

So far we have mainly focussed on hard sphere fluids 
near hard walls. In the following we discuss to which 
extent the aforementioned observations can be extended 
to other systems. This will be discussed along the lines 
of Fig. 5. 

Figure 6 displays the corresponding results for a fluid 
which is governed by a square-shoulder pair potential of 
finite strength (SUq, acting as a representative for inter¬ 
action potentials of finite range. The comparison with 
Fig. 5 does not reveal qualitative changes. Figures 7 and 
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(c) spherical wall 
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FIG. 9: Reduced coefficients 7n(<5), n £ {0,1, 2, 3}, characterizing the curvature expansion in Eq. (18) as function of the shift 
parameter J and as obtained from Eq. (25). Panels (a), (b), and (c) belong to the same systems as in Figs. 6(a), 7(a), and 8(a), 
respectively, with an additional soft part V of the external substrate potential: (a) V given by Eq. (Bl), (b) V given by Eq. (B2), 
and (c) V given by Eq. (B4). The data correspond to the following choice of system parameters: /3Vp'^ = pVp'^ = PVp'’ = 0.1 
(see Appendix B), Ly^/L = 1, and Ly-mlL = 5. 


8 , respectively, display the corresponding results for a 
fluid with a Yukawa pair potential, representing exponen¬ 
tially decaying interaction potentials, and for a Lennard- 
Jones pair potential, representing algebraically decaying 
pair potentials. Although comparing them with Figs. 5 
and 6 reveals certain differences, the main conclusions 
remain the same: the coefficients 7n(^) are strongly af¬ 
fected by the choice of the convention for <5 and, whereas 
for none of the systems considered here the curvature- 
dependence of the interfacial tension is exactly in agree¬ 
ment with morphometric thermodynamics, the morpho¬ 
metric form of the interfacial tension may be an excellent 
approximation for suitable conventions for 8. 

In order to further assess to which extent the above 
findings are generic, an additional excess part V of the 


external potential Y®’'* is considered. This excess part V 
is obtained by integrating the pair potential Uw{r) be¬ 
tween a fluid particle and a wall particle over the volume 
V of the wall: 

PV{v) = Qp, j d^r'pUp,{\v-v'\), (26) 

V 

where in V the number density Qw of the wall is taken 
to be constant. The particle-wall potentials Uu,{r) are 
chosen to be of a similar form as the pair potentials 
U{r) between the fluid particles (Eqs. (4)-(6)), with the 
exception that the Yukawa-like wall-particle potential 
is not truncated and the repulsive part of a Lennard- 
Jones-like wall potential is replaced by a hard wall. The 
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pair potentials are characterized by an energy scale Uqw, 
a length scale and a cut-off length L^w such that 
^ Lew) — O’ 

• the square-well (C/oiu < 0) or square-shoulder 
{Uqw > 0) potential 

PUUr < = pUow (27) 

with = Lew, 

• the Yukawa potential {Uqw > 0, Lew = oo) 

PUwir) = /317ou;— exp , (28) 

r \ LwJ 

• the Lennard-Jones potential (Uqw > 0) 

= -(%^)) ^Lw<r<Lew 

[oo ,r < Lw- 

(29) 


The resulting excess parts V of the substrate potentials 
as function of the radial distance z to a spherical reference 
surface of radius R are given in Eqs. (Bl), (B2), and (B4) 
in Appendix B. 

In Fig. 9 the results for systems with a non-vanishing 
excess part V of the external potential are shown. The 
parameters are chosen such that, apart from V ^ 0, the 
same systems as in Figs. 6 (a), 7(a), and 8 (a) are ana¬ 
lyzed. It again turns out that the coefficients 7 n(^) ex¬ 
hibit a strong dependence on S so that the quality of 
morphometric thermodynamics as an approximation de¬ 
pends sensitively on the convention for S. As a general 
trend, for all three examples shown in Fig. 9 the non- 
morphometric coefficient \15^^{5)\L'^ is not negligible in 
wider ranges of conventions for 5 than in the cases with¬ 
out an excess substrate potential V. In this sense the 
quality of morphometric thermodynamics as an approx¬ 
imation deteriorates in the presence of excess parts V of 
the substrate potential. 


via a square-well or square-shoulder potential (Fq. (4)), 
a Yukawa potential (Eq. (5)), or a Lennard-Jones poten¬ 
tial (Eq. ( 6 )). As shown in Figs. 5-8 the coefficients 7 n(^) 
of the curvature expansion in Eq. (18) may depend sen¬ 
sitively on the chosen interface convention, which is ex¬ 
pressed in terms of the shift parameter <5 (Fig. 1). There 
are conventions for which the morphometrically allowed 
coefficients are much larger than the morphometrically 
forbidden ones so that within them morphometric ther¬ 
modynamics is a reliable approximation of the interfacial 
tension. However, the opposite situation can occur for 
other interface conventions, in which case morphometric 
thermodynamics has to be used with caution. In par¬ 
ticular the reliability of morphometric thermodynamics 
as an approximation deteriorates in the presence of ex¬ 
cess contributions to the wall potential (Fig. 9). Based 
on these results, it turns out to be necessary in future 
applications of morphometric thermodynamics to clearly 
state which interface convention is chosen and why mor¬ 
phometric thermodynamics is expected to be a reliable 
approximation for that particular interface convention as 
compared with others. 


Appendix A: Ideal Gas 


In this appendix we analyze the exactly solvable case 
of non-interacting particles. The density functional in 
Eq. (7) with t/(r) =0 is minimized by the equilibrium 
number density 


eeq(r) = ^.^q-'^^exp (-^F'^"*(r)) with 

0X^ = A-^exp(/3fl). 

For pointlike ideal gas particles the convention J = 0 is 
convenient and will be used throughout this appendix. 
For this choice the interface of area A, the reference sur¬ 
face, and the geometrical wall surface are the same and 
the interfacial tension 7 is given by 


IV. CONCLUSIONS AND SUMMARY 

By using density functional theory within the second 
virial approximation we have analyzed several model flu¬ 
ids with small number densities in contact with planar, 
spherical, or cylindrical walls. We have compared the 
curvature expansion in Eq. (18) of the interfacial ten¬ 
sion 7 with the expression derived within morphometric 
thermodynamics (Eq. (3)). Particular attention has been 
paid to the implications of the choice of the position of the 
interface, which underlies the definition of the interfacial 
tension (Eq. (16) and Fig. 1). For none of the considered 
systems the expression for the interfacial tension in accor¬ 
dance with morphometric thermodynamics is exact, re¬ 
gardless of whether the particles interact with each other 


„bulk p 

h = —d\(exp(-^U(r)) - 1),(5 = 0 . (A2) 

The integration volume V in Eq. (A2) equals the volume 
accessible to the fluid particles. Therefore the integrand 
depends only on the excess part V of the external poten¬ 
tial U®’'* (see Eqs. (10) and (11)). 

In the case of a hard wall with V = 0 the interfacial 
tension of the ideal gas is zero, 7 = 0 , irrespective of 
the shape of the wall. Therefore the ideal gas is a useful 
choice for studying the influence of the excess part V ^ 0 
of an external potential on the morphometric coefficients. 
Here we analyze a Yukawa-like interaction Uw (Eq. (28)) 
between the fluid particles and the wall particles. The ex¬ 
cess part V follows from Eq. (26). For planar, spherical. 
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and cylindrical walls, respectively, one finds 


Viz) 

vr 

V{r) 

vr 


= exp - 


La, 


plane. 


=-exp - 


(A3) 


X [{S — l)e'^ + (S' + l)e '®] , sphere, 


V{r) 

vr 


=2Sh{S)Ko 


La, 


cylinder, 


where Vr = ‘^t^QwUowL^ denotes the strength of the ex¬ 
cess part V at contact with a planar wall and S = R/L^i 
with Kq and Ii as familiar modified Bessel functions. 
The expressions in Eq. (A3) for the excess parts V of the 
external potential facilitate to determine exactly the co¬ 
efficients jn of the curvature expansion of the interfacial 
tension 7 in Eq. (A2). In the case of a spherical wall we 
obtain 


r = 


bulk T 
0eQ 


00 

-E 

n—l 


{-pvrr 


nl n 


Lw 


- 1 — n 

n 


2 r ^2 


i ?2 


n" n 3 2 

r:::—I" TT ~ 1-1- n 


-^(2-n)(l-n) 


+0 


7^4 

zA. 

i?4 


O ( exp ( —2 


R 

La,, 


n n 
1 n 

"Eh 

,(5 = 0 , 

(A4) 


and for the cylindrical wall the corresponding result is 
given by 


— ^bulk T 

P7 — Peq ^ 


{-pvrY' 


i?2 


n—l 

1 1 

2n 8 


Lw 

'1 

I 

n 

R 

n 

2 

2 


+ O|§)},i = o. 

(A5) 


The expression for the planar wall is included in the ex¬ 
pressions for the curved walls (Eqs. (A4) and (A5)) as 


J 


the term being independent of the radius R. 

In Eqs. (A4) and (A5) the respective curvature expan¬ 
sions are presented up to and including the leading non- 
morphometric coefficients (belonging to R~^ in the case 
of spherical walls and to R~^ in the case of cylindrical 
walls) which in general are non-zero. 

Further interesting insight can be gained by studying ra¬ 
tios of particular coefficients: 


• The ratio of the leading coefficients 71 (belonging 
to R~Yj 



(A 6 ) 


i.e., a constant value independent of the strength 
PVr of the external potential. 

• For small amplitudes pvr the ratio of the subdom¬ 
inant coefficients 72 (belonging to i?“^) is given by 



(A7) 


For /3Vr ^ 1 this ratio reduces to the constant 
value 8/3 which has also been found in Ref. [11] 
and in the numerical calculations of Sec. III. 


• When comparing the leading ( 7 )^) and the subdomi¬ 
nant ( 7 ^) coefficients for cylindrical walls we obtain 

^ = -l^vr +^ ipvrf+o ({pvrf) . (A8) 

This implies that for <C 1 one has I 72 I 

1 7 /1 which contradicts the morphometric prediction 
according to which 72 should be zero. 

These results for an ideal gas of non-interacting parti¬ 
cles invalidate morphometric thermodynamics if there is 
a non-vanishing excess part ^ 0 of the external 

potential. 


Appendix B: Excess parts of the external potentials 


The calculation of the excess parts V of the external potentials according to Eq. (26) results in lengthy expressions. 
Here these are presented for spherical walls as function of the distance z > 0 from the reference surface of radius R, 
i.e., the distance from the center of the spherical wall is r = i? -|- z. 

In the case of the square-well or square-shoulder potential (Eq. (27)) the excess part V of the external potential is 
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given by 


pv{z > 0 ) = X < 


z ^ Lw • 
z<L^-2R: 

-2R < z < Lu 


/3^;« = -TTLipUo^g^. 


2 -^ 

Ll 

1 Li 
8 {R+zP 
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- + 6 ^ 




Lt 




-X?; (4f - 24£ + 8) + (-30f + 24£ - s) 


(Bl) 


-^(-12f+24f-6£) 


3 q 


In the case of the Yukawa potential (Eq. (28)) the excess part V is given by 


l3V{z > 0) =/?y/“ exp 

=27Tg^(]Uo^Ll 


Lyj ) R -\- z 


A 

Lon 


1 + exp 





(B2) 


(Note that Eq. (B2) and the sphere expression in Eq. (A3) are identical.) 

In the case of the Lennard-Jones-like potential (Eq. (29)) the reference surface and the geometrical wall surface do 
not coincide: R = Rg + L^, with Rg denoting the radius of the geometrical wall. Integration according to Eq. (26) 
results in an expression pVg{zg, Rg) in which Zg measures the distance from the geometrical wall surface, i.e., the 
distance from the center of the spherical wall is r = Rg + Zg\ 
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(B3) 


The excess part V of the external potential as function of the distance z from the reference surface of radius R is 
related to Ig in Eq. (B3) via 


l3V{z > 0) = pVgiz + L^,R- L^). 


(B4) 
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